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Abstract. Triggered di-hadron correlation studies using central Au+Au collisions at \J Snn 
= 200 GeV in STAR revealed a novel ridge-like structure in two dimensions (A??, A<f>) for 
high pr particles. A similar structure is also present in an inclusive un-triggered di-hadron 
correlation analysis. We study the (pr) evolution of di-hadron correlations by increasing the 
lower pr acceptance of both charged particles. A smooth evolution of data is observed and 
our results reproduce the triggered analysis structure near {pr} = 2.7 GeV/c. We quantify 
the correlation structure evolution by fitting a model function. The model function emphasizes 
possible initial state fluctuation contributions via the use of higher harmonic model components 
v n (n=l,2,3,4,5) and a remainder which is modeled via an asymmetric 2d Gaussian. The 
extracted parameters are compared to model predictions and p+p data at \/Snn = 200 GeV 
and possible origins of the nearside structure are discussed. 



1. Introduction 

Two dimensional, high px di-hadron correlation studies in STAR revealed a novel correlation 
structure elongated in pseudo-rapidity (Ar]) [1]. This structure which is known as the "ridge" 
was observed by applying an asymmetric transverse momentum cut on trigger (4 GeV/c > p? 

> 3 GeV/c) and associate (pt > 2GeV/c) particles [2]. Furthermore, it was revealed that the 
structure persists up to trigger particle pt > 6 GeV/c, possibly suggesting a jet-related origin. 
An independent analysis [3] also observed the ridge correlation for charged particles with px 

> 0.15 GeV/c for both trigger and associate particles. This analysis differs from the above 
mentioned one in three key aspects. It uses all possible correlated particle pairs, it applies a 
common lower pt threshold to all particles and the correlation function is normalized in a way 
that yields correlated particles per final state charged particle. The measured structure is not 
attributed to expected physical processes and in the past several attempts have been made to 
explain the long-range correlations through the medium modification of elementary processes 
such as jets or flux tube formation [4,5,6]. 

For our study we adopt the un-triggered analysis method and focus on the correlation structure 
in the 0-10% centrality bin where we have the highest probability of forming a deconfined 
medium of quarks and gluons. We then investigate the evolution of the correlation structure as 
a function of transverse momentum. The structure is studied by increasing the lower transverse 
momentum acceptance for both particles. The resulting correlation is modeled via an empirical 
model fit function. We discuss our findings in the context of models that describe medium 
properties and modified jet phenomena in deconfined QCD matter. 



2. Data and analysis 

The data used in this analysis were collected during two runs, in 2004 and 2009, using the STAR 
detector at the Relativistic Heavy Ion Collider (RHIC) Brookhaven National Lab (BNL), Long 
Island, New York. For the results presented here we have analyzed 11 million Au+Au central 
trigger events and 264 million p+p minimum bias events at \/Snn = 200 GeV. Only events with 
a primary vertex within |x,y| < 3cm and |z| < 25cm were taken into account. Charged particle 
tracks were reconstructed by the STAR TPC only and constrained to the primary vertex. The 
correlation function was calculated only for tracks within the kinematical cuts of |r/| < 1.0, full 
2tt acceptance in azimuth <fi and px = 0.15-15 GeV/c. The correlation measure for this analysis 
is defined as follows: 

Ap _ p s ib — Pref 
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The symbol p^ refers to the total correlated and uncorrelated pair density as a function 
of Arj and A<p, for hadrons from a single event. p re j refers to the number of uncorrelated 
pairs as a function of Ar] and A<f>, for hadrons from different events. Ap is therefore the 
number of correlated pairs in a single event. The denominator is the square root of the number 
of background pairs, which corresponds to the number of particles. The correlation function 
therefore measures the number of correlated pairs per particle. In order to correct for acceptance 
effects, we sub-bin the z vertex into ten 5cm bins and allow a multiplicity window of 50 for event 
mixing. The conversion electron/positron background is reduced by using a 1.5cr dE/dx cut on 
electrons in the momentum ranges, 0.2 < pr < 0.45 GeV/c and 0.7 < px < 0.8 GeV/c. In order 
to analyze the transverse momentum evolution of di-hadron correlations in the 0-10% centrality 
bin we steadily increase the lower particle pr to gain further insight into the origin of the long 
range correlations. The correlation function evolution for eight selected bins is shown in Fig. [TJ 
We model structures via the use of an empirical function and infer potential underlying physics 
mechanisms of the near side long range correlation. In the following section we discuss the model 
fit function. 

3. Model function study 

The main part of the model function is based on recent theory developments related to initial 
energy density fluctuation contributions to the two particle correlation structure [7] . It has been 
pointed out that the fluctuation contributions are not negligible even in the raw data spectrum 
and should be observable in very central events using high px charged particles, as shown in 
Fig. [2] The remaining structure after subtracting out the higher order Fourier harmonics, which 
we term the remainder, is modeled via an asymmetric 2d Gaussian, and thus can be attributed 
to modified jet phenomena. The remainder properties are quantified by comparing to p+p data. 

3.1. Evidence for initial density fluctuations in data 

Following the theory prediction [7], we focus our attention to the 0-1% centrality bin. Fig. [2] 
shows the projected 2d correlation structure on to A(f> as a function of pt revealing a double 
hump structure in the 2 < px < 5 GeV/c interval for the away side correlations. This observation 
is consistent with the kinematic ranges predicted by theory and most likely due to triangular 
flow [8]. 

3.2. Model function 

We initially attempt to describe the data using a maximum number of Fourier coefficients, which 
for our data requires coefficients up to n = 5 in the cos(nA^>) series. However the remaining 
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Figure 1. Evolution of the raw normalized correlation structure for selected cuts on px of the 
particles (increasing (pr) from left to right) at 0—10% centrality bin. 
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Figure 2. Angular correlations as a function of pt cuts on the pairs for 0-1% central data. 



structure reveals a significant residual, thus requiring an additional asymmetric 2d Gaussian 
term. The complete function can be broken down as follows. 



F = c oc Offset (h) 

+ ci * cos(lA0) oc Away-side momentum conservation (fi) 

+ C2* cos(2A(j)) oc Second order Fourier harmonic (fy) 

+ C3 * cos(3A(j)) oc Third order Fourier harmonic (f±) 

+ C4 * cos(4A(f>) oc Fourth order Fourier harmonic (fe) 

+ C5 * cos(5A0) oc Fifth order Fourier harmonic (f%) 
+ cq * exp(—0.5* ((A(j)/cj) 2 + (Ar]/c&) 2 )) oc Modified jet fragmentation (fj) 



F = /1 + /2 + /3 + /4 + /5 + /6 + /7 
In the following section we discuss and interpret the model parameters we extract. 

4. Results and Discussion 

4-1- Results 

In this paper we mainly focus on the evolution of the model components shown in Fig. |3j 



Fig. 3(a) shows the evolution of the remainder component after subtracting out the harmonics 
(n=l-5) contributions to the correlation structure. Fig. |3(b)| shows the sum of all harmonics 
(n=2-5) contributions. The remainder parameters might shed some light on potential modifi- 
cation to jets in the medium whereas the harmonics parameters can be used to test predictions 
from hydrodynamical models. 

The Ar]-A(f) asymmetry of the remainder widths are shown in Fig. |4} We observe an 
asymmetry below (px) ~ 2.5 GeV/c. The largest asymmetry is observed at pt ~ 1-3 GeV/c 
where the Ar] width is a factor of two greater than the Acfi width. 

We compare the extracted parameters to the same parameters from our analysis of p+p data 
in Fig. [5] In the region above 2.5 GeV/c where the Arj-A(f> widths are symmetric in Au+Au, 
the widths are about a factor of two larger than in p+p. Below that (px) value, we observe 
modification in both Ar] and A<f> remainder withs compared to corresponding p+p widths. 
Furthermore, the maximum relative width modulation (AA/pp) is a factor of two greater in Ar] 
width compared to A<j) width with the maximum modulations occuring at two different (pt) 
values. 

Fig. [6] shows the comparison between the remainder and the sum of all harmonics amplitudes. 
The amplitudes are comparable in the 0.9 < pt < 2.1 GeV/c range. The summed harmonics 
amplitude, which extracts the near-side long range correlation in our model, saturates at a 
small value at high px- The correlation strength of the added harmonic structure peaks at a 
{pt) ~ 0.9 GeV/c which is higher than the (px) for bulk particles thus indicating that long- 
range correlations might carry particles coming from jet fragmentation (hard processes). The 
asymmetric remainder amplitude drops with (px) and the symmetric remainder amplitude rises. 
The behavior of the symmetric remainder amplitude at px > 2.5 GeV/c can be understood based 
on more correlated jet fragment pairs in a jet like cone. 

Finally we show scaling relations using higher order harmonics in order to compare our 
findings to certain hydro predictions [9]. Plotted in fig. 7 are the ratios v l J A jv^ 2 and v\ jv^ 2 . 
Theory predicts those ratios should be independent of {px) in the bulk pj- regime. Even 
though we observe the predicted trends, in order to extract more information about the medium 
properties from the scaling ratios, we need to compare our results to theory applied to similar 
kinematic constraints (centrality and px) than the data. 



(a) The remainder evolution as a function of pr 
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(b) V2 + v 3 + f 4 + v$ evolution as a function of pr 



Figure 3. The evolution of the "remainder" (a) and the summed v n structure (n = 2 - 5) (b), 
as a function of the lower cut on px of the particles. 
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Figure 4. The (pr) evolution of Ar/, Acj) widths of the remainder. 
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Figure 5. Comparison of the remainder Ar], A(j) widths to corresponding p+p widths. 



4-2. Discussion 

We studied the evolution of di-hadron correlations in angular variables Ar/, Acj) using Au+Au 
collisions at \/Snn = 200 GeV. Our main focus was to study the long-range correlation observed 
in the central events bin by increasing the lower pj* acceptance of both charged particles. By 
focusing on 0-10% centrality bin, we study the most violent heavy ion collisions where a formation 
of a quark gluon plasma is expected. Therefore our study of the long-range correlation could 
reveal further insight into the properties of the deconfined medium. 
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Figure 6. Remainder amplitude variation compared to summed harmonic (v2 + V3 + V4 + U5) 
amplitude. 
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Figure 7. Higher order harmonic scaling relations using, ^2,^3 and U4. 



The measured structures evolve smoothly as a function of (pr) ■ The near-side correlation evolves 
from an elongated 2d Gaussian-like structure towards a structure which has a flat pedestal and a 
jet like peak on top of it. The away side structure broadens and flattens at higher px- Both the 
near side long-range correlation and away side amplitude reduces with increasing px- However, 
we observe a jet like peak strength growth (see fig. [T]) after (pr) 2.5 GeV/c possibly due to 
an increase of correlated pairs within a fragmenting high px parton even though the total yield 
in the correlation function (including uncorrelated background) drops, as expected. 
The data evolution was modeled using an empirical model function which has nine free 
parameters. We extract the contribution to the correlations due to initial energy density 
fluctuations via the use of higher order Fourier harmonics cos(nA0), n = 1,2,3,4,5. In addition, 
the data require an asymmetric 2d Gaussian model component in order to fit the residual 
structure after subtracting the harmonics contribution. 

The evidence for initial energy density fluctuations in the correlation data (see Fig. implies 



that we do create a medium that behaves collectively. We further test the collectivity of the 
medium by testing our predictions against hydrodynamical scaling relations (see Fig. [7]) and 
observe the predicted trends. Theory states that any deviation of the scaling ratios from ideal 
hydro predictions could be due to partial thermalization effects [9]. However our present scaling 
data could not be directly compared to those predictions due to discrepancies in kinematic cuts. 
Nonetheless, our extracted v n parameters are useful in studying medium properties(e.g. viscosity 
and diffusion coefficients) and initial conditions (e.g. CGC vs. Glauber initial distributions) of 
the colliding heavy ions. 

The remainder component after subtracting harmonics, can be viewed in two parts. Below (pr) 
~ 2.5 GeV/c the remainder widths are asymmetric and the amplitude drops as a function of 
(jpr)- I n the same (px) region, the widths are modified compared to p+p data. The relative 
asymmetries as well as width modulation with respect to p+p data occur at different (pr) values. 
Nonetheless, the remainder below (pr) ~ 2.5 GeV/c is elongated in At/ and more modified 
than the Acj) width in comparison to p+p parameters. Therefore the cause of the asymmetric 
remainder component is open to interpretation and the measured modulation factors could be 
useful input for jet shape modification models. On the other hand, the symmetric remainder 
component above {pr) ~ 2.5 GeV/c shows a width modification of a factor two compared to 
p+p collisions. Theorists could use this information to constrain partonic energy loss models 
and calculation of energy transport coefficients in the deconfined medium. 
In summary, we discussed that the near side long range correlation structure can be understood 
as arising from initial energy density fluctuations and possible modified jet phenomena. The 
modified asymmetric 2d Gaussian needs to be studied further in the context of jet fragmentation. 
Charge dependent and particle identified measurements are two suggested approaches for further 
studies. 
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